arXiv:1509.03376vl [physics.optics] 11 Sep 2015 


Inhomogeneity-related cutoff dependence of the Casimir energy and stress 


F. Bao,i’2 j. S. Evans,2 M. Fang,^ and S. 

^Department of Physics, Zhejiang University, Hangzhou 310058, China 
^Centre for Optical and Electromagnetic Research, 

JORCEP, Zhejiang University, Hangzhou 310058, China 
^ Centre for Optical and Electromagnetic Research, ZJU-SCNU Joint Research Center of Photonics, 
South China Normal University, Cuangzhou 510000, China 
'^Department of Electromagnetic Engineering, Royal Institute of Technology, lOOjf Stockholm, Sweden 

(Dated: September 14, 2015) 

The cutoff dependence of the Casimir energy and stress is studied using the Green’s function 
method for a system that is piecewise-smoothly inhomogeneous along one dimension. The asymp¬ 
totic cylinder kernel expansions of the energy and stress are obtained, with some extra cutoff terms 
that are induced by the inhomogeneity. Introducing interfaces to the system one by one shows how 
those cutoff terms emerge and illuminates their physical interpretations. Based on that, we propose 
a subtraction scheme to address the problem of the remaining cutoff dependence in the Casimir 
stress in an inhomogeneous medium, and show that the nontouching Casimir force between two 
separated bodies is cutoff independent. The cancellation of the electric and magnetic contributions 
to the surface divergence near a perfectly conducting wall is found to be incomplete in the case of 
inhomogeneity. 


I. INTRODUCTION 

Decades ago Casimir theoretically predicted the attrac¬ 
tion between two neutral, perfectly conducting plates by 
virtue of the zero-point energy [1]. His seminal work 
started the exploration of a new territory [2-4] where 
the vacuum fluctuations result in previously unexpected 
effects. The Casimir force was confirmed by later exper¬ 
iments [5-8] and might find applications in nano/micro- 
electromechanical-systems (NEMS/MEMS) [9-12]. De¬ 
spite the wide interest in the Casimir effect, due to its 
fascinating physical effects such as the quantum repulsion 
[13-15] and levitation [16-19], quantum friction [20, 21], 
quantum torque [22, 23], and also due to the deep rela¬ 
tionship between vacuum energy and the cosmology con¬ 
stant [24], there remains open questions regarding the 
cutoff dependence or the divergence in the Casimir en¬ 
ergy and stress. 

The cutoff terms which need to be renormalized trace 
back to Casimir’s original work where he considers the to¬ 
tal vacuum energy and subtracts a term that is already 
present even without the plates. This term is not subject 
to the boundary conditions and turns out to be diverg¬ 
ing. A similar treatment can be found in the general 
theory [25] of van der Waals force (also Casimir force) by 
Dzyaloshinskii et al., where they have subtracted from 
the stress a term that one would get if the medium is un¬ 
bounded, uniform and of the same local property e. This 
subtraction is based on the assumption that the short 
waves do not sense the “inhomogeneity” (changing of e 
across different objects) and thus do not contribute to the 
Casimir force. Casimir and Dzyaloshinskii’s treatments 
have successfully predicted a force that agrees quite well 
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with experiments, and thus are widely adopted as an es¬ 
sential procedure (subtraction of the “empty space”, the 
“unbound medium”, or the “bare/bulk contribution”). 
However, this is not always adequate. When the space 
is abnormal [14], or curved geometries are present [26- 
28], the energy or stress will still be diverging after the 
essential procedure of subtraction, even with the disper¬ 
sion considered [29]. The remaining cutoff dependences 
make it ambiguous to calculate some physical quantities 
(the value of the cutoff needs to be determined by exper¬ 
iments) like the surface tension associated with the sur¬ 
face divergence [30], and some of the cutoff terms even 
do not seem renormalizable [27, 31-33]. Fortunately, the 
Casimir force between separated bodies is not influenced 
and is always finite (free of cutoff dependence), see the 
discussions in [34] and in Sec. 4.3.3 in Ref. [4]. This has 
led to many successful methods [35-40] which concern 
merely the force between separated bodies enclosed by 
arbitrary manifolds. 

Recently, the remaining divergence in the stress ten¬ 
sor has been reported in the case where two flat half 
spaces are separated by a continuously inhomogeneous 
medium [41, 42], even with dispersion considered. Cu¬ 
riosities mainly lie in the question whether the Casimir 
force between separated bodies is cutoff dependent or not 
and the problem on how to obtain a finite and meaningful 
stress to describe the Casimir interaction. Nevertheless, 
allowing for inhomogeneity of the medium not only com¬ 
plicates the nontouching Casimir interaction, but also 
permits the nonvanishing local pressure from the medium 
on a single compact body, and moreover, implies the non¬ 
uniformity of the surface tension at the interfaces. The 
electromagnetic contribution to the local pressure and 
the surface tension is cutoff dependent, which is reason¬ 
able as it actually relates to the Casimir self interaction 
and the interaction between two touching objects—the 
compact body and the surrounding medium. The usual 
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law of the Casimir force a~^ implies this cutoff depen¬ 
dence when the separation a —?► 0. In a previous pa¬ 
per [43], by considering an infinitesimally thin film so as 
to neglect the cutoff-dependent pressure [44] , we have 
shown that the nontouching Casimir force (NCF) is cut¬ 
off independent, in the first-order perturbation theory. In 
this paper, we ignore the dispersion of the media, adopt 
the Green’s function method to study the cylinder kernel 
expansions of the energy and stress, and propose a modi¬ 
fied Dzyaloshinskii’s subtraction [25] to retrieve the non¬ 
touching Casimir interaction (note that the NCF is not 
an independent observable unless the nonuniform pres¬ 
sure balances out and the surface energy is conserved). 
The cutoff dependences are analysed term by term and 
the unresolved renormalization of the quadratic and log¬ 
arithmic divergences found in curved geometries are dis¬ 
cussed here in the context of inhomogeneity. 

This paper is organised as follows. In Sec. II, we show 
in detail how we can get the asymptotic expansion of 
the Casimir energy, for different setups. In Sec. Ill, we 
discuss how the cutoff terms change when the setup is 
modified so that we can assign the physical meanings 
reasonably. A modified Dzyaloshinskii’s subtraction is 
also given. In Sec. IV, we study the cutoff dependence of 
the Casimir energy and stress for two systems of different 
kinds of optical response. The cancellation between the 
electric and magnetic contributions to surface divergence 
near a perfectly conducting (PEC) wall is discussed. We 
also show numerically that the stress is finite based on 
our modified subtraction scheme. And we conclude in 
Sec. V. Throughout this paper, we use the natural units 
H = c = 1. 


II. METHODS 

We study the fluctuating electromagnetic field in a 
medium that is inhomogeneous along one dimension, say 
the X axis. The stress tensor and the energy density can 
be rewritten in terms of the Green’s tensor. Due to the 
translation invariance along the y-z plane, the Green’s 
tensor can be reduced to two scalar Green’s functions 
which account for the electric and magnetic contribu¬ 
tions, respectively [41] 


importance of the expansion terms, and then give the full 
results of the electromagnetic problem in Sec. IV. 

The scalar field equation is [9^ — V • e“^V] ip = 0, and 
its Hamiltonian is 

n^^dt(pdt(p+^e~^^(f^(f. (2) 

For a scalar field interacting with a background by the 
Lagrangian Cint = — we refer the readers to 

Ref. [34]. With point-splitting regularization, the vac¬ 
uum or Casimir energy density is 

^=( 01 -^ 10 ) 

= hni i [dtd' + e-^V-V] {0\Tip(t,vMt\r')\0) 

2 (3) 

= lim ^ [dtd'^ + • V'] A(t, r; t', r'), 

where T is the time-ordering symbol and A is the Feyn¬ 
man propagator that satisfies [i9( — V • V] A = S{t — 
t')S{r — r'). For stationary boundaries, we can use the 
time-translation invariance to obtain the Fourier image of 
A, G(r, r',aj), satisfying [— — V • e“^V] G = S(r — r'). 
Performing Wick’s rotation u = ik, finally we have 

OO 

J [-K^ + e-iy • V'] GdK. (4) 
0 

Hereafter, the limit r' — r is implied in the expression 
of the Casimir energy density. For a dielectric function 
that varies along the x axis, using the space-translation 
invariance along y-z plane, the energy density could be 
simplified further as 

OO 

U = jj ue~^ -\-u^ dx&^ ( 5 ) 

0 

where Q{x,x' ,u,iK) meets exactly the bottom definition 
in Eq. (1) with p = 1. 

Consider the differential equation 


[^£K^-f ^Ge) =S{x - x'), 

[USK^ -I- - da:e~^dxe]{e~^GM) =S(x - x), 

where y{x) is the permeability, e{x) is the permittivity, 
u is the wave number along the y-z plane, dx is the full 
derivative with respect to x, and k is the imaginary fre¬ 
quency w = iK. In this section, we give a general mathe¬ 
matical description of these two scalar Green’s functions, 
and seek the asymptotic expansions of the Casimir en¬ 
ergy density and stress. For simplicity, we start with a 
scalar version of the problem, corresponding to the mag¬ 
netic contribution when /i = I, to highlight the physical 


[e - dxe-^dxe] V = 0, (6) 

where = v? -\- eK^,£ > 0 and u = kcos6. When 
e is smooth and = y/u^ >> 1 thus the 

wave number fc ^ 1 (hereafter, fc ^ I always means 
» 1), Y can be approximated by the WKB method 

to arbitrary order, Y = exp | i ^”S'„(u, k, a:)|. 

Substituting into the equation yields two types of 
V, one of which is A = fi{u, and 

the other is i? = fi{u, k, where /i = 

exp {lEfci, 3.5... .5**5'*} and /2 = | 0 . 2 . 4 .... 



3 


Specifically, we have 


^0 = ^k, 

Si = -^{2 + sm^ 9), 

4:6 


S 2 S 0 = 


SaSi = 


r(12-5siii^6»)- 


4 (2 — sin^ 0) e 


SaSI = 


32e2^ ' 23 e’ 

24+12sin2 6»-15sin®6»£3 
26 £3 

16 + 4 sin^ 0 — 9 sin'* 9 ee 
25 ? 

2 — sin^ 9 £^3) 

1008 + 960 sin^ 9 + 600 sin'* 0-1105 sin® 0 ¥ 


248 + 196 sii+ 0 + 50 sin^ 0 - 221 sin® 0 £^£ 
2® £3 

lO + 5sin2 0-7sin'*0££(3) 


25 


in £2 


28 + 12 sin 0—19 sin 0 e 
27 

2 - sin^ 0 £("*) 

25 T’ 


( 7 ) 


which will be used later. Here dot means the full deriva¬ 
tive with respect to x. The exact form of these two so¬ 
lutions is usually unavailable and thus is less interesting 
to us. Instead, we aim to construct the Green’s function 
Q using them. The result is 

Q = Qbl + Gb2 + Qs, 

= ( 0 (x-x')-i), 

< . 1 ( 8 ) 

Qb2 = -^[AB' + J!B]. 

^ . 

Here W = AB — Ab = — 2 / 4/2 is the Wronskian that has 
the property dx {eW) = 0, 0 is the Heaviside step func¬ 
tion, Cl and C 2 are arbitrary constants which are deter¬ 
mined by the boundary conditions, A = eA, A' = e'A, 
and so on [in the following contents, A = dx{eA)\. Gs 
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is sometimes referred to as the scattering Green’s func¬ 
tion, and Gb = Gbi + Gb 2 referred to as the bare Green’s 
function. Although Gbi alone is a special solution to the 
point-source equation that G obeys, it is Gb that is the 
well-defined singular part of G, because Gbi and Gb 2 both 
diverges when k^oo or x — x' ^00 but Gb does not. 
When the medium is homogeneous, say vacuum, we have 
A = ■^e^^,B = _ _ 2 ^ ^] 4 gj 4 have 

Gb = ^ which is exactly the vacuum Green’s 

function. When the medium is smoothly inhomogeneous, 
Gb differs from the bare Green’s function for the homoge¬ 
neous medium of the same local e [45] , and consequently 
the energy density also differs. 


A. Smooth inhomogeneity 


If the inhomogeneity we consider is smooth and the 
system is unbounded, Ci and C 2 are both zero. When 
substituted into Eq. (5), Gbi always gives a contribution 
of S{x — x') to the kernel of the integral. The energy 
density reads 



As can be seen from the WKB approximation, when 
/c » 1, we can expand AB as ^ = 

YAi=T o,i{x,9)k~’‘. Specifically, we have = 0 for 
all i € Evens, ai = 03 = ^SoS2, and 05 = 

I — SqS^J ■ Also we have 

. . . . ( 10 ) 
(SoS2y-2S^S3S + 2S^S4 , ^„_ 4 , 
p k 

where S = Si+dxlne. Introducing an upper cutoff of the 
wave number A = S(x—x') to the integral (note the lower 
bound is not zero because we require k :$> I), or alterna¬ 
tively using the cutoff parameter t = I/A = min{x —x'}, 
and leaving aside the cumbersome calculations, we have 
the regularized Casimir energy density 


,reg _ 1 ,-4 , 1 504££ - 509£^ n 

6720£^ 

1 II52£3£G) _ 5280£2££(3) _ 4240£2£2 + 15680££2£ - 7155£* 

~ 47r2^ 30720£4 


Int + 0 ( 1 ). 


( 11 ) 


The leading term is the energy density for a homoge¬ 
neous medium, which can be easily verified by the mode- 


summation method 

,7 


1 /■*/* dk3 1 

¥7eJo (^ 16 ^^ 


( 12 ) 
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The subleading terms are nonzero when the medium is 
inhomogeneous. 


B. One interface 

Consider the inhomogeneity which is smooth every¬ 
where except at X = 0. Denote L the left domain x < 0, 
I the interface, and R the right domain x > 0. When 
X, x' > 0, along with the direct propagation from x' to 
X, now the scattering takes place = Gb + Gf- Using 
the continuity of G and e~^dxG at the interface, we can 
obtain from Eq. (8) the scattering Green’s function 


a) from the right (RL) or left (LR) side. In the other 
domain where x, x' < 0, the scattering Green’s function 
is just like the above with all L and R interchanged, and 
all A and B interchanged. For clarity, the superscripts in 
the following formulas will be suppressed if they can be 
figured out by the context. 

The quantity r is actually the generalised reflection 
coefficient at the interface. Expanded when fc ^ 1, 
gives a leading term that is identical to the p-polarised 
Fresnel formula 


r(o) = _ 


e^k^ +e^k^ 


(15) 




where 


r 


RL 

0 


e^dx — e^dx In^-^ 

e^dx — e^dx 


(13) 


(14) 


Hereafter subscript 0 (or a in the following contents) 
means the quantity is evaluated in the limit x —>■ 0 (or 


Its first subleading term is 


2e^k^{e^S^ - 
{e^k^ + 


(16) 


The second and third subleading terms are given in 
Appendix A. All the higher terms are unnecessary for 
analysing the cutoff dependence. 

Substituting the scattering Green’s function into 
Eq. (5) and letting x —>■ x', we get 





^ B ^BB{Ao AB 


(17) 


where 


AB Bo 


k^ - 2kS -f 52 -h 2 S 0 S 2 


^{SoS2S + S^S3) + 0{k-^) 


g-2/2|S^ 


(18) 


and = e xhe exponent can be expanded as 

- 2 / 2 IS = -2r'^o|5 - 2^521? - 2^^541? + 0(e'^). (19) 

The leading term —2^“^5'o|g = —2 k{x)dx will be kept 
in the exponent while the subleading terms can be pulled 
down by an exponential expansion. Recall the expansion 
of ^ and rg, we see there are four different types of k 
that have been used in the energy density—fcoi k{x), k{x) 
and To unify the notation, we define 


/ si x') 

k{x) = 1 /cos^ 9o -\ - - sin^ Oq x ko = J\f{x, 0)ko, 

y £0 

' kix) = \ cos^ 9o + sin^ 9o x ko = Mix, 0)ko, 

V eo 

= \/cos^ 9o + — siM 9o X ko = Mioo, 0)ko, 

v V £0 

( 20 ) 

where M indicates the transformation of the wave num¬ 
ber (an analogue of the refractive index). Now the energy 


density becomes 

5 /*7r/2 poo 

dkoF{9o,x)d9o + 0(1), 

n=l -^0 

( 21 ) 

where T is some function of 9o and x, and X = 
A/’(x,0)dx is an analogue of the optical distance. 
When X —0 so A —>• 0, the k integral in Eq. (21) 
yields five diverging terms X~^, X~^ and In A 

(precisely these five terms come with a factor 
but this factor doesn’t affect the diverging property; the 
last term In A is integrable near the interface). The 
divergence of the energy density when the interface is 
approached, the so-called surface divergence, is a well- 
known phenomenon [26, 46, 47]. Specifically, the term 
A“^ relates to while the term A“^ involves up to rU) 
and thus . Therefore, the surface divergence would 
disappear if the continuity of at the interface is 

met. A similar result has been reported by K. Milton in 
Ref. [48] . There, the scalar field in a monomial potential 
[d^ — -|- x“]G = i5 is studied. And the surface diver- 





















5 


gence is found to disappear when a > 2 (though the two 
results are slightly different because of the different mod¬ 
els. They consistently show that, the surface divergence 
will disappear if the medium is smooth enough). 

It is still challenging to get an explicit full expansion of 
the energy density for an arbitrary permittivity. The dif¬ 
ficulty lies in getting function J^. To simplify the problem 
while still focusing on the influence of the inhomogeneity, 
we keep only the first two terms in the above exponential 
expansion, so that exp{— 2 / 2 |g} = — 2 ^S' 2 |g). 

By doing this, we are able to get the first subleading 
term of the expansion of the energy density. Besides, 
we let £'^ —?► oo so that the reflection coefficient is also 
simplified 


1 -I- 0dx^^e)QX, X = a;[l -|- |(sin^ lne)oa;] and 

<S'2|o = {S 2 )ox. At last, the energy density is 


ui 


1 

Idn^y/eo 





If it is oo, we have ^ O S and L ^ R those 

interchanges in Eq. (17). But in an elegant way we will 
have the same result as Eq. (23). 


C. Two interfaces 


= -1 - 25o/fco. (22) 

Furthermore, after the k integral in Eq. (21), we take 
the near-interface limit a; 1 so that Af{x, 0) = 


When the medium contains two interfaces at x = 0, a, 
letting L, C and R denote the left, central and right do¬ 
mains respectively, the scattering Green’s function in the 
central domain takes the form 


rLR 

ys 


1 1 
eW 1 - 


AA' 



+ BB' 



(AS'-f A'S)iPr“rf^| , 


(24) 


where V = (^)q describes the one-round 

propagation of the photon in the central domain. 

Physically, the photon propagator from source x' to 
point X where 0 < x,x' < a, includes the contributions 
of direct propagation and reflections by the interfaces. 
The former gives the bare Green’s function, while the 
latter gives the scattering Green’s function. However, 
the single reflection by one interface accounts for the in¬ 
teraction between two adjoining domains, and does not 


contribute to the nontouching interaction between two 
separated objects which intuitively should involve mul¬ 
tiple reflections between two interfaces. Therefore, from 
Eq. (24) we should subtract Eq. (13) and a counterpart 
for the interface at a, i.e., Gc = Gg^ — Gg— Gf- Then we 
get the Green’s function that describes the nontouching 
Gasimir interaction we are interested in 


Gc 


1 1 
eWl- 


AA' 



The factor V is essential to make the energy density Uc 
(yielded by Gc) finite near the interface. To see this, we 
don’t need to do the asymptotic expansion, but note that 
it is proportional to exp [—2 fcdx], an exponentially 
decaying term that makes any polynomial integrable over 
k. It follows that, the nontouching Gasimir force is cutoff 
independent. 


III. ANALYSES 

By using the macroscopic quantity e to describe the 
medium and its interaction with the field, we actually 
imply the possibility of self interaction of the medium 
mediated by the field. A well-known picture of that is, 




the medium can be polarized by a photon by an operator 
n = — r'), and the free energy is thought of as 

the interacting energy between polaritons mediated by 
virtual photons [25]. Not only the dielectrics, even the 
vacuum can produce electron-positron pairs and is thus 
polarizable. In light of these, the Gasimir energy given by 
Eq. (3) or (29) can be understood as a part of self energy 
of the medium. Specifically, according to the previous 
section, when the medium is smooth and unbounded, the 
energy density is Ub- Thus it corresponds to the self in¬ 
teraction of an unbounded smooth medium. When there 
is one interface, the energy density is = Ub+uf'^^'^. 

The newly arising energy thus corresponds to the 

interaction of two adjoining domains. When there are 


- {AB' + A!B) 




(25) 
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two interfaces, the energy density is 

= Ub+U^ +U^ +Uc. (26) 

The energy lAc corresponds to the nontouching interac¬ 
tion between two separated domains. 

Let us focus on the case of two interfaces. The first 
three terms of Eq. (26) are all cutoff dependent. To 
remove the infinities, Dzyaloshinskii et al. proposed to 
subtract a term that we would get if the medium is un¬ 
bounded and homogeneous with the same local property 
e, that is, the leading term in Eq. (11), based on the as¬ 
sumption that short waves do not sense the inhomogene¬ 
ity and thus do not contribute to the Casimir force. This 
is reasonable for real materials, and at least works for ho¬ 
mogeneous media. In the latter case, all the subleading 
terms in Eq. (11) vanish, and is independent of the 

position thus vanishes in calculating the force, leaving us 
only Uc that we are actually interested in [Note that if we 
are considering the EM field in a homogeneous medium, 
Ug -I- (and similarly stress) will vanish near a PEC 
boundary due to the cancellation between the electric and 
magnetic contributions, as shown in Eq. (32)]. Again, in 
the case of the homogeneous medium, the Dzyaloshin- 
skii’s subtraction has removed all the energy that is not 
subject to the boundary conditions (the bare Green’s 
function contribution), equivalent with Casimir’s origi¬ 
nal treatment. 

When the medium is smoothly inhomogeneous and 
non-dispersive, arbitrarily short waves contribute to the 
Casimir force, so the Dzyaloshinskii’s subtraction fails. 
In order to fix this, Philbin et al. proposed removing the 
contributions of the small-scale inhomogeneity by using 
the Green’s function given by the first-order WKB ap¬ 
proximation as Eqs. (36) and (37) in Ref. [41] (up to So 
and Si). Here we see, this cannot remove completely all 
cutoff terms caused by Qb for an arbitrary e, but indeed 
relaxes the divergences. To remove Ub from U, we need 
higher-order approximations up to iS' 4 . 

The last term of Eq. (26) is called the free Casimir 
energy, because it can be released by changing the sepa¬ 
ration distance, giving rise to the nontouching Casimir 
force. However, when the medium is inhomogeneous, 
we see the energy will also change along with the 

shift of the interface, invalidating the usual definition of 
the free Casimir energy as the energy difference when 
the separation is finite and when the separation is in¬ 
finite Ec = E{a) — E{oo) (invalidating Casimir’s treat¬ 
ment). Consequently, the nontouching Casimir force is 
no longer an independent observable unless the integral 
of the scattering energy over all interfaces, here the sum 
U^+U^, gives a constant. Notably, for the EM field, the 
electric and the magnetic contributions now do not can¬ 
cel completely for PEC boundaries, in contrast to the 
case of a homogeneous medium. To retrieve the free 
Casimir energy (stress), we propose the following sub¬ 
traction scheme 

(27) 


or one can directly use the corresponding Green’s func¬ 
tion as Eq. (25). 

When the energy density is integrated over the whole 
space, we get the well-known cylinder kernel expansion 
of the Casimir energy [33]. Homogeneous (sectionally) 
medium and flat boundaries would yield the simplest re¬ 
sult E ^ Vt~'^ + St~^ + Ec + 0{t) where the inverse 
cubic term can still be zero when some kind of cancella¬ 
tion takes place. The presence of curvature of the bound¬ 
aries [26, 28, 29] would immediately introduce some ex¬ 
tra terms so that E ~ Vt~^ + St~^ + e 2 t~^ + 63 ^“^ -I- 
f\nt + Ec + 0{t). Here all the coefficients are related 
to the geometric characteristics of the system. Particu¬ 
larly, V is the volume of the space, S is the surface area 
of the boundaries, and the others have more complex 
meanings. According to dimensional analysis, 62 should 
be some length, 63 is a pure number, and / is the in¬ 
verse of a length. To renormalize these cutoff dependent 
energies, one could introduce the classical energy of the 
system = pV + ^S + Ee 2 + ke 3 + hf and absorb t in 

these classical quantities p, 7 , F, k and h [4, 27]. However, 
an unresolved problem is, there seems to be no standard 
physical meanings for F, k and h (p being the pressure 
and 7 being the surface tension). 

Here in our case, the interfaces are flat, and these ex¬ 
tra cutoff terms are caused by the inhomogeneity of the 
medium. Integrating the self energy Ub over the space, 
we won’t have an explicit factor V which we would have if 
the medium is homogeneous (and we can have if we use 
some averaged quantity), and also there are quadratic 
and logarithmic terms. Besides, the integration of the 
interacting energy Ug would yield S{t~^ +t~^ + t~^ + In t) 
but not only St~^. In this way, we see that, though the 
volume term Vt~'^ must still come from Ub and the sur¬ 
face term St~^ must be still from Ug, the other terms 
especially the quadratic and logarithmic terms are also 
from Ub and Ug , and contain contributions from both of 
them. Note that those terms involve higher-order deriva¬ 
tives of £, thus together with S or under the volume inte¬ 
gration, they have coefficients of indeed the dimensions 
same as 62 , 63 or /, respectively. 

Based on the above analyses, we want to take a 
heuristic interpretation of these cutoff terms—terms in 
62 ^”^ + eot~^ -f /Int are not to be renormalized by some 
new classical quantities, but are only higher-order correc¬ 
tions to the pressure p and the surface tension 7 . That 
is, introduce = f^pdV + fg^'jdS, and absorb the 

cutoff terms in Ub by p and in Ug (precisely f Ugdx) by 
7 . The physics is clear. The self interaction of one object 
contributes to its inner pressure, and the interaction of 
two adjoining objects contributes to the surface tension 
at the interface. 

At last, we should mention that, the energy density 
exactly at the interface is not well-defined, though the 
Green’s function is still continuous there. This is ob¬ 
vious as can be seen in Eq. (11) —to define the energy 
density we need the continuity of up to . This might 
be overcome by softening the interface layer by specify- 


Uc=U^^ EUb-U^ -U^, 
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ing a smooth profile of the permittivity. It needs future 
study. 

IV. APPLICATIONS 

In this section we apply the above method to the elec¬ 
tromagnetic field and give the asymptotic expansions of 
the energy density and stress, to discuss some inter¬ 
esting topics. Notably, the Minkowski’s stress tensor 


= DiEj+BiHj — ^5ij (D • E -k B • H) and the Raabe- 
Welsch stress [50] , the ordinary energy density of the elec¬ 
tromagnetic field U = \ (D • E -f B • H)^ d^ and the 
Casimir energy density [51] derived from the interaction 
of dipoles all show similar asymptotic behaviours with 
minor differences in the coefficients. Therefore, we give 
only the results of the Minkowski’s stress and the energy 
density derived from it. The component of the stress that 
generates force is 


oo 7r/2 

^—= J k'^dk J ^^l~^[k'^ - d^d'JGE + £~^[k‘^ - cos 


47r2^ 


eae, 


and the energy density is 


0 0 


oo 'K ]2 


u = 


47r2 0Ie , 


k^dk / ^[k^{cos^ 9 — siTL 9)dxd'^\GE 


0 0 


-|-£’ 


"^[fc^(cos^ 9 — sin^ 9) + cos0d0. 


(28) 


(29) 


Those two scalar Green’s functions are given in Eq. (1) 
and have similar solutions as Eq. (8) with /i or £ in ap¬ 
propriate places. Note also that for electric and mag¬ 
netic Green’s functions we have different A's and B's. 
In principle, arbitrary permittivity and permeability that 
vary along one dimension are tractable, but the results 
are sufficiently cumbersome that only by adding further 
constrains can we get relatively simple analytical expres¬ 
sions. Two examples we will discuss below are the trivial 
magnetic response system fj, = 1 and the constant light 
speed system fie = 1. In the results, we keep only the 
first two terms of the expansions. 


A. Trivial magnetic response system 


When /r = I, the WKB approximation for the mag¬ 
netic Green’s function is the same as Eq. (7), while for 
the electric Green’s function we have 


50 = ^k, 

51 = sin^ 9, 

4e 


(30) 


The bare parts of the regularized energy density and 


stress are 


?/reg _ _ 


t-^ + 


1 224££ - 229£^ 




req 

=- 


=t-4- 


3360£2 
23£2 


t ^ -k O(lnt), 


247r2^ 47r2^480£2 


t 2 _|_ O(lllt). 


(31) 


When there is one interface at a: = 0 and £l —>■ oo 
in this or the constant light speed system, the electric 
Green’s function obeys the Dirichlet boundary condition 
while the magnetic Green’s function obeys the Neumann 
condition. The latter condition yields the same reflec¬ 
tion coefficient as Eq. (22), while the Dirichlet one yields 
exactly te = 4. The fact that te has an opposite sign 
against tm, to some extent, results in the cancellation 
between the electric and magnetic contributions. The 
scattering parts of energy and stress read 


ui 


a 


R 

S 


-I 

I 

Idn'^y/e^ 




(32) 


The absence of the term in the energy density is 
due to the above mentioned EM cancellation, while we 
see when the medium is inhomogeneous, the cancellation 
is incomplete. The absence of that term in the stress, 
however, is a generic property, because the operator dxd'^ 
acting on the scattering Green’s function would produce 
a leading term so that to the leading term, (/c^ — dxd'^) 
















FIG. 1. The dotted red curve recovers the previous result [41]. 
While using our proposed subtraction scheme, we obtain the 
solid green curve. The dashed blue curve is the stress when 
£ = 2. The inset is the profile of e. 


Making a substitution e ^ e we can get the approxi¬ 
mation of the electric Green’s function expressed by e 


r = 1, 


*■ = 5- 

+1- 


(34) 


The bare parts of the regularized energy density and 
stress are 


1 1 £2 

7/res _ I ^_ 

'' 87r2 ^ 47r2 24^2 


t ^ + O(lnt), 

2 (35) 




The scattering parts are 


vanishes in the expression of the stress. Again, we see 
when the medium is inhomogeneous, there are subleading 
terms. And, this is exactly the reason to the divergence 
near the interfaces reported in [41]. As analysed in the 
previous section, here we give a numeric illustration of 
the proposed subtraction scheme [using Eqs. (25) and 
(28)] to fix the problem raised by the above reference. 

As can be seen in Fig. 1, the profile of the permittivity 
is plotted in the inset [e^ = = 1, and = 3e~^ 

(solid) or = 2 (dashed)]. The dotted red curve in 
the main figure is simulated with — Qb (the previous 
method, where B and A are the first and second kind 
of the modified Bessel functions, respectively) which still 
contains -I- so that the stress diverges near the 
interfaces. Using C/c, immediately we arrive at the solid 
green curve. The dashed blue curve is for the case = 2. 
Note that, either of these two satisfies the Casimir 
repulsion criterion thus it’s reasonable to 

see the solid green curve in the same side (negative side 
of the vertical axis) with the dashed blue curve (while the 
dotted red curve is in the positive side that corresponds 
to attraction). 


B. Constant light speed system 

When = 1, = 1. Thus for the magnetic Green’s 

function, we have 


r = 1 , 


=-^ 


^ 2^0 = 


2 £’ 

3e^ £ 

8e2 Ae' 


(33) 


ui 


a 


R 

s 


-1 

167r2 

1 

167r2 




(36) 


This constant light speed system is interesting be¬ 
cause, the usual surface divergence near a curved bound¬ 
ary was reported to disappear [52, 53]. Here, we see 
the inhomogeneity-induced surface divergences are still 
present in this constant light speed system. 


V. CONCLUSIONS 

To summarise, we have studied the asymptotic cylin¬ 
der kernel expansion of the Casimir energy and stress 
within an smoothly inhomogeneous medium, using the 
WKB approximation of the Green’s function. When the 
medium consists of several parts and thus is piecewise- 
smoothly inhomogeneous, we have interpreted the self en¬ 
ergy of the medium as possessed by the interaction among 
those parts. A heuristic renormalization of the cutoff 
terms has been given. Also, we have given a modified 
subtraction scheme to retrieve a finite stress to describe 
the nontouching Casimir interaction which consequently 
is finite as well, with numeric illustration. Moreover, we 
have studied the cancellation between the electric and 
magnetic contributions to the surface divergence near a 
PEC wall, and the surface divergence for the constant 
light speed system. For the former, we have found the 
cancellation is incomplete when the medium is inhomo¬ 
geneous. And for the latter, we have found the surface 
divergence is still present, in contrast to the case of cur¬ 
vature. 
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Appendix A: Reflection coefficients 


Following Eq. (16), the second order correction of the reflection coefficient is 


( 2 ) _ - e^S^)l _ .r 


{e^k^ + £^k^) 


(Al) 


0 


^0 


and the third order reads 

2£^k^{£^S^ - £^S^)l 24(e-^5^ - e^5^)o - S^S^) - S§{e^S^ + e^S^) 


^(3) ^ 


(e^fc^ + £^k^)Q 

O^LuRtoRqL qR ^RqLQLqR\ 

zepfeoi^e jg Jq Jq t jg Jq Jq Jo 
l,Li.R(^RUL I pLhR\2 

n.g n.Q r\j C rt 


{e^k^ + £^/c^)q 


(A2) 
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